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Abstract 



We consider the classical linear stability of a static universe filled with a non-interacting 
mixture of isotropic radiation and a ghost scalar field. Unlike the conventional Einstein static 
O , model, this cosmology is stable against homogeneous and isotropic perturbations. This is 

shown by means of exact oscillatory solutions about the original static state. We also examine 



5-H , the linear response of the static ghost universe to all types of inhomogeneous fluctuations, 

namely density, vorticity and gravitational-wave perturbations. The results show that the 
static background constrains the linear evolution of these distortions, to the extent that 
' density perturbations remain time invariant, vortical distortions vanish and gravitational 

^ , waves oscillate with constant amplitude. We discuss the potential implications of these 

' results for past-eternal initial states in classical general relativistic cosmology. 

m 

• ■ 1 Introduction 
O 

I Cosmologists have retained an interest in the consequences of admitting 'ghost' scalar fields into 

C5 ' the Universe. These fields may have a negative energy density [I], and are far from new, having 

been incorporated in the early study of steady-state cosmology O [3], in the investigation of 
scalar fields [4], phantom dark matter OE], and K-essence [7]. Although ghost-fields may be 
' unstable at the quantum level [8l[9], they provide a simple theoretical laboratory for exploring 

the physical consequences of cyclic closed universes which bounce at a finite radius. In par- 
ticular, Barrow et al [TO] examined the unusual consequences that arise for theories of varying 
'constants' in such oscillating universes. Ellis and collaborators IT^ have also considered 
cosmologies that emerge into expansion from a past eternal Einstein static state, reminiscent of 
the original Eddington-Lemaitre cosmology favoured by Eddington as an infinitely old universe 
that can nonetheless be thermodynamically young \13\ I14j. One further consequence of these 
ghost cosmologies is the existence of static solutions with new properties. In this paper we will 
examine the stability of Einstein static universes in these theories. Two new features will be 
of interest. First, the existence of almost static solutions that oscillate about a finite radius; 
and second, the classical linear stability of certain ghost-like static models under all three types 
of inhomogeneous perturbations (i.e. scalar, vector and tensor modes). The inhomogeneous 
situation shares some analogies with the stability shown by the Einstein static universe in the 
presence of conventional (non-ghost) matter fields, as discussed in [15j-[18jP] 



^The stability analysis discussed there is also rendered delicate because the Einstein static universe has compact 
space sections and Killing vectors and so is a conical point in the space of all solutions to the Einstein equations 
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Throughout this paper we assume a pseudo-Riemannian spacetime with a Lorentzian metric 
of signature (— H — We also adopt the 1+3 covariant approach to general relativity and 
cosmology, referring the reader to older and more recent reviews of the formalism for further 
reading and details [!22]-|24j. 

2 Ghost scalar fields 

At the basis of the 1+3 covariant formalism is the concept of the fundamental observers, which 
when introduced allow for a unique 'threading' of the spacetime into time and 3-dimensional 
space. The worldlines of these observers are tangent to a timelike 4-velocity field, relative to 
which physical quantities and equations decompose into their timelike and spacelike components. 
The formalism combines mathematical compactness with physical transparency and has been 
used to a variety of cosmological are general relativistic studies. Here we will use the covariant 
approach to study cosmological models containing a ghost scalar field (tp) with zero potential 
(i.e. V{^) = 0) and a Lagrangian of the form [lOj 

>C^ = ^V,VV>V=5, (1) 

where g is the determinant of the spacetime metric (gab) and is the standard covariant 
derivative operator. Then, the associated stress-energy tensor reads 

^if = - Va^VfeV' + I Vci^V'i'9ab ■ (2) 

To achieve an 1+3 decomposition of the ^-field, we assume that VaV' is timelike (namely 
that V(j'0V"'i/' < 0) over a given spacetime region. In that case is normal to the spacelike 
hypersurfaces ipi^a) = constant and defines our fundamental timelike 4-velocity field 

Ua = --^Vall^, (3) 

with ip = u'^Vaip 7^ 0- This guarantees that tp'^ = — VaV'V"'0 > and that Uau"' = — 1 
as required. The above 4-velocity also defines the fundamental time direction and the tensor 
hab = dab + UaUfy that projects into the 3-space orthogonal to Ua- The projection tensor also 
provides the covariant derivative = ha^Vhj which operates in the observers' instantaneous 
rest space. 

In covariant terms, the kinematics of the Ua-congruence ^ are determined by means of the 
irreducible decomposition [I3l [24l 

Vfetia = ^ Qhab + (Tab + ^^ab " ^a«6 , (4) 

where 0, aab and uJab respectively describe the average volume expansion/contraction, shear 
distortions and the vorticity of the fiow. Following our 4-velocity choice, the later vanishes 
identically (i.e. uiab = 0), making the Ua-field irrotationalU The 4-acceleration vector {Aa) 

which exhibit linearisation instability about it unless higher-order constraints are imposed |19)-|21). 

^An direct consequence of our 4-velocity choice (see ^) is that Daip ~ 0. This in turn guarantees that 
i^ab = T)[tUa] = and the irrotational nature of the flow (e.g. see [24]). 
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describes non-gravitational/non-inertial forces and is given by 

Aa = -iBai;. (5) 

The above, which may be seen as the momentum conservation law of the ip-field (compare to 
Eq. (j31b ) in ^ 14. 2p . imphes that ip acts as an acceleration potential. 

Introducing the Ua-frame ([3]) also facilitates a convenient fluid-like description of the ip- 
field. In particular, with respect to the fundamental observers, the energy-momentum tensor 
([2]) recasts into 

T^T =P^'^^UaUb+p^'^hab, (6) 

where 

pW=pW = -i^^ (7) 

The right-hand side of ([7|) is negative by default, ensuring that p^^^ < and consequently the 
ghost nature of the '0-field. 

In the presence of a non-vanishing potential, scalar fields do not generally behave like 
barotropic fluids. Here, however, the potential has been set to zero and a barotropic description 
of the V-field is therefore possible. According to ([7]), our scalar field corresponds to a 'stiff' 
medium with effective barotropic index 

/p(V') = I and a corresponding sound-speed 

given by ci = dp^"^^ /dp^"^^ = 1. 



3 Homogeneous ghost cosmologies 
3.1 Static spatially closed models 

The evolution of a Friedmann-Robertson- Walker (FRW) cosmology with nonzero spatial curva- 
ture and a non-vanishing cosmological constant is monitored by the set 

H' = lip + A)-^, H = -H'-^(p + 3p) + h (8) 

and 

p = -3H{p + p). (9) 

In the above p and p represent the total energy density and pressure of the matter, K = 0, ±1 
is the 3-curvature index and A is the cosmological constant. In addition H = a/a \s the Hubble 
parameter of the model and a is the cosmological scale factor. To close this system one needs 
to supply an (effective) equation of state for the (total) matter sources. 

Consider a static cosmology with zero cosmological constant and positive spatial curvature. 
Then, the cosmological scale factor is time-independent {a = a^ = constant) and can be seen as 
the radius of the universe. Setting A = 0, X = +1 and H = = H , Eqs. ([8l)-(l9]) recast into the 
constrains 

3 

/9=-^, p + 3p = and P = , (10) 

«o 
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respectively^ The former of the above shows that the size of our model is inversely proportional 
to the energy density of its matter content, namely that oq = \/3/p. The second guarantees 
that the effective barotropic index of the system is w = p/p = —1/3. 

Suppose now that matter forms non-interacting mixture of isotropic radiation, with p^"'^ = 
p^^^ /3, supplemented by a ghost scalar field (-0) with zero potential and an effective equation 
of state given by ([7]). When both components are comoving, their total energy density and 
pressure, relative to a fundamental observer, are p = p^'^'^ + p^^'^ and p = p^'^'^ -\-pW respectively. 
In that case expressions ([TO]) reduce to 



p(7)+pW=3^ p(7)+2pW = (11) 

and 

p(y) + = . (12) 

Of the above constraints, pTh ) combines with (flTb ) to ensure that the energy densities of the 
two species are directly related, according to 

pW = ApM = _p_ (13) 
This result then, together with (fTTb ). leads to 



W = ^^ao = y'-^ and P''' = ^ = J ^ (l^) 



expressing the radius of our static model in terms of the energy densities of the individual matter 
components. We finally note that, for non-interacting fluids, the conservation law (jl2p is also 
separately satisfied (i.e. p^'^^ = = p^^'^). 



3.2 Almost-static oscillatory models 

Before looking into inhomogeneous cosmologies with ghost-like scalar fields, we will provide an 
example of a soluble almost-static homogeneous model. The Friedmann equation for the scale 
factor, a{t), of a closed universe containing a ghost scalar field with density p^^^ = = 
—^a~^ < 0, and radiation with density p^"'^ = Ta~^, reads 

_ ^ T 1 

where ^ and T are positive constants. Expressed in terms of conformal time (??, with dry = 
a^^dt), the above integrates to give 

a^irj) = ^ {r + Vr2 - 4^' sin [2(7? + %)]} , (16) 

with r]Q constant when > 4^. 

■^A positive 3-curvature allows us to relate more directly with the recent Emergent-universe scenario [111 112) . 
as well as the much earlier Einstein static and the Eddington-Lemaitre models [l3]-[T8]. A nonzero cosmological 
constant is not necessary for our purposes, since its role is now played by the dynamical, ghost-like scalar field. 
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Identifying the global maximum and the global minimum of the expansion with a^j^^ — 
[r + -4^']/2 and a^j^ = [r - Vr^ -4^']/2 respectively, we find that the scale factor 
evolution is given by 

a^iv) = \ {oLx + Omin + (aLx - amin) sm [2(r/ + 7?o)]} . (17) 

Given that both Omax and Omin are finite (recall that T, ^ ^ 0), this solution oscillates about an 
also finite radius between successive maxima and minimao. If we had chosen = then there 
would have been a single minimum (and no maximum) with collapse for a < a^im and expansion 
for a > amin- Since tp = \/2^a~^, the scalar field evolution given by [10] 



tp = ±2tan"^ 

With Cmax ^ amin; we have amin = \/^/r and Omax = VT- We can then see that the bounce 
duration is At ~ a^jj^/amax- Since tjj ~ \/6r^/2/^ near the bounce, we find that Aip ~ \/6, 
independently of initial conditions, during each bounce. 

The extreme case = 4\I' is a static universe with = a^ax — ^min = r/2 = V^. Setting 
d = and a = 0, we can see that this case is realized when = —pr/2, giving a = y/6/pr- 
We can see explicitly that the solution = r/2 is stable against homogeneous and isotropic 
conformal perturbations that do not change the curvature of the model. Indeed, deviations from 
the static model are monitored by the parameter 

e = VT^ - 4^ = al^, - a^, . (19) 

Then, following (jl6p and (jl7p . the perturbed spacetime oscillates about the static solution with 
(arbitrary) amplitude determined by e. In particular, expressed in terms of the above defined 
parameter, solutions p^ . (fT7|) assume the perturbative form 

«'(^?) = ^{r + esin[2(r? + r?o)]} , (20) 

which describes an almost static oscillatory (or cyclic) universe of constant amplitude. This 
result is in sharp contrast to the instability shown by the conventional Einstein-static cosmol- 

ogy m-M- 

We will now turn to consider the linear stability of these static solutions against all three 
types of inhomogeneous perturbations namely scalar, vector and pure-tensor (i.e. gravitational- 
wave) distortions. In the latter case we will also take a step into the nonlinear regime, by looking 
at the static model's response to gravity-wave anisotropics of arbitrary magnidude. 

4 Perturbed static ghost cosmologies 
4.1 Describing inhomogeneity 

In accord with the 1+3 covariant approach to cosmology, spatial inhomogeneities in the density 
distribution of the matter fields are described via the associated comoving fractional density 

*For mathematical simplicity we have specialised the matter content of the model to be radiation plus the 
ghost field. The qualitative behavior of the Friedmann equation would be similar if the radiation was replaced by 
any perfect fluid with p > p > — p/3. 



r tan (?? + r]o) + VT^ - 4^' 
2\/^ 



(18) 
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gradients. For the total fluid, these are given by the dimensionless ratio |25j 



a 



Aa = -BaP, (21) 
P 

with Da = ha^'Vi, representing the covariant derivative operator in the observers rest space. 
Similarly, the gradients 

Ai'^) = -^D,pW and A^^^ = ^bJ^\ (22) 

monitor inhomogeneities in the densities of the two individual matter components. When the 
background model is spatially homogeneous all of the above gradients vanish and therefore 
describe linear density perturbations in a gauge invariant manner ^26j . Also, definitions (j21|) 
and ([22]) imply that 

A. = i(,(-)AW+pWAW). (23) 

The density gradients are supplemented by additional inhomogeneity variables, of which the 
most important describes fluctuations in the volume expansion (0). The latter is represented 
by the orthogonally projected gradient 

Za = aD,e , (24) 

which is also independent of the gauge choice. 

Consider the background model described in § 13.11 and perturb it by introducing weak inho- 
mogeneities and anisotropies. To linear order, definitions ([2T|) . (p2]) and (p^ read 

«0 „ A (-y) ^0 J-y) A (^) «0 



A - ^ D o A(^) - — ^ D o(^) A^'^) - — ^ D o^^^ (25) 

and 



Za = aoBaQ , (26) 
respectively. Also to first order, relation (f23l) reduces to 

A, = 2AW - Ai'^) , (27) 

since p^'^'^ = —p^"'^/2 = —p in the background (see conditions pT]) ). 

4.2 Linear density inhomogeneities 

In non-interacting multi-component systems, the individual members observe separate conser- 
vation laws. Since the V'- field does not interact with the radiative component, their associated 
linear momentum conservation formulae are 

ao{l + w^^)Aa = -cf^A^\ (28) 

with i = 'J, ip and Aa being the common 4-acceleration vector (e.g. see [21])- In the absence of 
interactions, linear density inhomogeneities in each one of the species evolve according to 

A« = -(l + u;«)z,, (29) 
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where 

= - i p (l + Scf) A, - ^ p (1 + 3u;)^a + aoBaB'^Ab , (30) 

to linear order. Note that p, p, w and correspond to the total fluid, and is given by (j2ip . 

When dealing with a perturbed static universe filled with a non-interacting mixture of ra- 
diation and a ghost scalar field, we have w^'^^ = 1/3 = (?s^'^^ and w^'^-* = 1 = Cs^'^^ at the zero 
perturbative level. Hence, the linear expressions ([29]) - ([28]) lead to 



AaoAa = -AW , 2aoAa = -A^ (31) 

and 

Ai^) = -^^a, AW = -2Z„. (32) 

The first of these sets implies that aI''^^ = 2Ai'^\ while second ensures that A^'^^ = 2Ai'^'*/3 + C 
(with C being time- independent). Therefore, to linear order 

aW=2AW = ^C. (33) 

In other words, linear density inhomogeneities (both in the individual species and the total fluid 
- see Eq. (|27p ) remain constant in time. This result seems to reflect the high symmetry of 
the unperturbed model and the properties of the two matter flelds. In particular, the static 
background, the absence of a cosmological constant and the non-interacting nature of the ma- 
terial mixture mean that the energy densities of the constituent species are directly related (see 
constraints (fTT|) - (fT^ ). When the two matter flelds also share the same 4- velocity, as it happens 
here, the aforementioned constraints lead to the momentum conservation laws (j31|) FI The latter, 
together with the propagation formulae ([32p. guarantee that any linear density inhomogeneities 
that might exist will not change in timeO 

Result ([33|) also guarantees the time independence of linear density perturbations (described 
via the scalar A = aoD^'Aa) as well as of shape distortions (monitored by means of the trace-free 
tensor A^^^^ = aoDip/S.^)). On the other hand, vortical perturbation in the density distribution 
of the species (described via the antisymmetric tensor Aj^^ = aoDj^Aa]) vanish identically. To 
verify this recall that (e.g. see Appendix A. 3 in |24j ) 

D[6Da]P = pwab = 0, (34) 

with the null result guaranteed on all perturbative levels by the irrotational nature of the Ua- 
frame (see § [2] earlier) . 



4.3 Linear gravitational-wave perturbations 

So far, our analysis has shown that scalar (density) perturbations are time independent, while 
vector (rotational) distortions vanish. The remaining third type, concerns pure tensor fluctua- 
tions. Gravitational waves are covariantly described by the transverse parts of the electric {Eab) 

^In general, one could assign different 4-velocities to the two fluids, in which case relative-motion terms will 
also appear in the linear equations (e.g. see § 2.4 and § 3.3 in [24]'). 

®The time independence of Ai''' and aI*' implies that the linear expansion inhomogeneities must vanish at 
all times (see Eqs. (132} ). This imposes a strict condition to the right-hand side of (|30[) . which must remain zero 
at all times. Checking the consistency of the linear condition goes beyond the scope of this article. 
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and the magnetic (Hah) components of the Weyl tensor. To hnear order, these are isolated by 
switching off both the scalar and the vector perturbations (to guarantee that T)''Eab = = T)^Hab 
- e.g. see [21] )■ On a static, FRW-type background, the linear propagation of these two fields is 
monitored by the set 

Eab = --^ P(l + w)aab + curlif^b (35) 

and 

Hab = -cuvlEab ■ (36) 

where curli^afe = £cd{a^'^Eb)'^ for any traceless and transverse tensor Eab- Also to first order, the 
magnetic part of the Weyl field is related to the shear by means of 

3 

CnilHab = ^CTab- D^CJafe , (37) 

with = D^Da representing the orthogonally projected covariant Laplacian and ag the radius 
of the background 3-D hypersurfaces. The above ensures the direct interconnection between the 
magnetic Weyl field and the shear at the linear perturbative level. This in turn allows us to 
replace Eq. (j36p with the shear propagation formula 

O'ab = —Eab ■ (38) 
The latter combines with (|35p to give the wave-equation 

1 3 

0"a6 = /0(1 + w)aab 2 ^"-b + D^C^afe , (39) 

Z Oq 

of the transverse part of the shear tensor. When the static background contains a mixture of 
non- interacting radiation and a ghost scalar field, w = —1/3 and p = S/oq (see § 13. ip . Applied 
to this environment, Eq. (|39p reduces to 

2 

f^ab = 2 + D^CTai, . (40) 

Our next step is to Fourier decompose the last expression. Introducing the standard ten- 

(n) (n) 

sor harmonic functions, we may write aab = J2n'^in)Qab ^ where Q^^^ are the aforementioned 
harmonics with Q^^ = qJ^^^, si^^ = = D'-Qi^ ancQ 

D^Ql? = -(^)'Qlt (41) 

Substituting decomposition aab = Z]n<^(n)Qib^ i^ito Eq. (00]), the harmonics decouple and we 
arrive at the wave equation of the n-th shear mode 

^(") = ~ ~ ^{n) ' (42) 



where > 3 due to the compactness of the unperturbed spatial sections. Expression (j42p 
governs the linear evolution of gravitational waves on a static, FRW-type background that has 



^Angled brackets denote the symmetric and trace-free part of orthogonally projected second rank tensors. 
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spherical spatial geometry and contains a mixture of non-interacting black-body radiation with 
a ghost scalar field. The solution of (j42p has the simple oscillatory form 



o"(n) = Ci sin 



( 



\/n2-F2 
al 



) 



t +C2 cos 



( 



«o 



) 



t 



(43) 



which ensures that gravitational-wave perturbation oscillate with a time-independent amplitude. 
In other words, the pure tensor modes are neutrally stable at the linear level on all scales. 

5 Non-linear gravitational-wave anisotropies 

The addition of large amplitude gravitational wave anisotropies can render the static universe 
unstable. This can be seen by examining Eq. (|15p for the isotropic and homogeneous background. 
Suppose that we add a simple anisotropy to the expansion, with isotropic 3-curvature. This is 
of the Bianchi I type and adds a new anisotropy energy density term for the shear, = S^/a^ 
with S constant, to the Friedmann equation, where a{t) is now the geometric mean scale factor 
and the Friedmann equation becomes 



Clearly, if the anisotropy is sufficiently large, so that > then neither the static solution nor 
the exact solution with finite non-singular oscillations around it ()16p - (jl7p any longer exist. This 
situation requires anisotropic perturbations to the FRW background that are far from linear. 

We can envisage an irregular early state in which the ghost field, radiation, anisotropy and 3- 
curvature index, K, all vary in space. The effective Friedmann equation for the spatially-varying 
scale factor a{x, t) would have the form 



Regions where S^(x) > ^'(x) would not oscillate through successive non-singular cycles; regions 
where k{x) < would bounce only once and then expand forever. Only those regions with 
low anisotropy and positive curvature would evolve through a succession of periodic cycles. 
If there was production of radiation entropy from cycle to cycle, in accord with the Second 
Law of thermodynamics, then this would be equivalent to an increase in the value of F from 
cycle to cycle [27\ and a corresponding increase in the amplitude of the oscillations around the 
static universe in Eq. p6|) . since they scale with F. They would approach the zero curvature 
expansion from below. If a positive cosmological constant were added to the right-hand side 
of Eq. (|45|) . then, no matter how small its magnitude, it would eventually come to control 
the dynamics if successive cycles grow in amplitude. Thereafter, the oscillations would cease 
and the expansion would asymptote towards an ever-expanding spatially-flat de Sitter state, as 
shown in [27]. More complicated forms of anisotropy, including 3-curvature anisotropies, would 
display a similar instability because at early times they include also the simple anisotropic stress 
considered here. 




(44) 




(45) 
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6 Discussion 



We have shown that in simple closed Friedmannian cosmologies containing black-body radiation 
and a ghost scalar field, there exist almost static solutions, which are oscillating about a finite 
radius. This is in sharp contrast to the instability shown by the Einstein-static cosmology, when 
conventional matter is used [13j . Our static universe also appears stable against all three type 
of linear inhomogeneous perturbations, in a way analogous to that seen in [l5]-[T8]. However, 
if large-amplitude homogeneous and anisotropic perturbations are introduced, the model be- 
comes unstable (it no longer exists), with its past and future evolution moving towards a strong 
curvature singularity. In the absence of these large perturbations the static solutions with sta- 
ble bounded oscillations could provide a viable past-eternal state for the Universe, although it 
would be subject to evolution if entropy was required to increase with time. Our analysis is 
of interest for studies of the very early universe incorporating theories of particle physics with 
new types of ghost matter field. It shows that it is theoretically possible for universes to get 
trapped in early-time evolutionary tracks that exhibit bounded oscillations about a non-singular 
static state. Similar behaviour might also arise in higher-order gravity theories [2H]-[2n], which 
introduce stresses that mimic the presence of ghost fields and would prevent some domains from 
participating in inflation. 
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